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We present measurements of the branching fractions for the charmless two-body decays B° — ► 
7r + 7r~ and B° — > K + -k~ , and a search for the decay B° — ► K + K~ . We include the effects of final- 
state radiation from the daughter mesons for the first time, and quote branching fractions for the 
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inclusive processes B° h+ti- 717, where h and h are pions or kaons. The maximum value of the 
sum of the energies of the n undetected photons, _E™ ax , is mode-dependent. Using a data sample 
of approximately 227 million Y(AS) — > BB decays collected with the BABAR detector at the PEP-II 
asymmetric-energy e + e~ collider at SLAC, we measure: 

150 MeV) = (5.1 ±0.4 ±0.2) x 10" 6 , 
105 MeV) = (18.1 ±0.6 ±0.6) x 10" 6 , 
= 59 MeV) < 0.5 x 10~ 6 (90% confidence level), 

where the first uncertainty is statistical and the second is systematic. Theoretical calculations can 
be used to extrapolate from the above measurements the non-radiative branching fractions, B°. 
Using one such calculation, we find: 

B°(B° -> TT+O = (5.5 ±0.4 ±0.3) x 10~ 6 , 
B°{B° -» K + tt~) = (19.1 ± 0.6 ± 0.6) x 10" 6 , 
B°(B° -> K + K~) < 0.5 x 10~ 6 (90% confidence level). 

Meaningful comparison between theory and experiment, as well as combination of measurements 
from different experiments, can be performed only in terms of these non-radiative quantities. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30.Er 



B(B° -» 7T+7V- n T , £™ ax = 
B{B° -f K + 7T~ n T , E™* = 
B(B° -> K+K- n T , E™** = 



Charmless hadronic two-body B decays to pions and 
kaons provide a wealth of information on CP violation 
in the B system, including all angles of the unitarity tri- 
angle. The time-dependent CP asymmetries in the tttt 
system can be used to estimate the angle a [l[ ; the decay 
rates for the Kir channels provide information on 7 [2|. 
Recently, direct CP violation in decay was established in 
the B system through observation of a significant rate 
asymmetry between B° — > K + tt~ and B° — > K~tt + de- 
cays [H, HI • As B physics experiments accumulate much 
larger data sets, charmless two-body B decays will con- 
tinue to play a fundamental role in testing the stan- 
dard model description of CP violation. Measurements 
of branching fractions for all the charmless two-body de- 
cays are invaluable in testing the various theoretical ap- 
proaches to the underlying hadron dynamics [f|- We 
present measurements of branching fractions for the de- 
cays B° — > 7r + 7r~ and A' + 7r~ Q, and a search for the 
decay B° — > AT + A'~ using a data sample about 2.5 times 
larger than that used for the most precise, previously 
published measurements @, H, [§] of these quantities. 

As radiative corrections have already proved to be im- 
portant in precise determinations of interesting quantities 
in the context of kaon physics [Io| , we account for them 
in this analysis as well. We can relate the observable 
decay rates T hh > (A" lax ) for B° — > h + h'~n"/ (and thus 
the branching fractions) to the theoretical non-radiative 
widths ry, ,, using the energy-dependent correction fac- 
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tors G hh ,{E™*^) 

= T hh ,(n)G hh ,(E™;ri, (1) 

where A™ ax is the maximum value allowed for the sum of 
the undetected photon energies and (i is the renormaliza- 
tion scale at which r* , and G hh > (£™ ax ) are calculated 
(the product being independent of fi). Extracting 
allows a more meaningful comparison with theoretical 
calculations and also between different experimental re- 
sults. Additionally, for A" lax at the kinematic limit, G 
approaches unity (to order ckqed / 7r )j so that the r° /; 
and the corresponding branching fractions, can be inter- 
preted theoretically in a cleaner way. 

The data sample used for this analysis contains 
(226.6 ± 2.5) x 10 6 TUS) -> BB decays collected by 
the BABAR detector at the SLAC PEP-II e+e~ 
asymmetric-energy storage ring. The primary detector 
components used in the analysis are a charged-particle 
tracking system consisting of a five-layer silicon vertex 
detector and a 40-layer drift chamber surrounded by a 
1.5-T solenoidal magnet, an electromagnetic calorime- 
ter comprising 6580 CsI(Tl) crystals, and a dedicated 
particle-identification system consisting of a detector of 
internally reflected Cherenkov light providing at least 3 a 
K-tt separation over the range of laboratory momentum 
relevant for this study (1.5-4.5 GeV/c). 

The data sample used in this analysis is similar to that 
used in the BABAR measurements of direct CP violation in 
B° — > K + tt~ Q and time-dependent CP-violating asym- 
metry amplitudes S^n and C^ in B° — > tt + tt~ [13j (the 
reader is referred to those references for further details of 
the analysis technique). Event selection criteria are iden- 
tical to those used in the CP analyses 0, EH , except that 
we remove the requirement on the difference in the decay 
times (At) between the two B mesons in order to mini- 
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Mode 


MC 


QED calculation 


7T + 7T~ 


89.8 ±0.1 


88.8 ±0.5 


K+TY- 


92.4 ±0.1 


91.7 ±0.5 


K+K~ 


94.7 ±0.1 


94.7 ±0.5 



TABLE I: Percentage of events with |AJ5| < 150 MeV and 
photon energy below the cut-off (2.6 MeV) in the Monte Carlo 
simulation, as given by the (see Tab.|I]) simulation and by the 
QED calculation described in the text. 



mize systematic uncertainties on the branching fraction 
measurements. 

We identify B° — > h + h'~ (h, hi — it or K) candidates 
with selection requirements on track and Cherenkov an- 
gle (0 C ) quality, B decay kinematic variables, and event 
topology. The final sample contains 69264 events and 
is defined by requirements on two kinematic variables: 
(1) the difference AE = E* B — \/s/2 between the recon- 
structed energy of the B candidate in the e + e~ center- 
of-mass (CM) frame and \fs/2; and (2) the beam-energy 
substituted mass toes = \/ (s/2 ± pi ■ pu) 2 /Ef — Pb 2 - 
Here, y/s is the total CM energy, and the B momen- 
tum pb and the four-momentum (E^pi) of the e + e~ 
initial state are defined in the laboratory frame. To sim- 
plify the analysis, we use the pion mass for all tracks 
in the track reconstruction and the calculation of the 
kinematic variables. We select those B candidates with 
\AE\ < 150 MeV, and 5.20 < m ES < 5.29GeV/c 2 . 

The efficiencies of the selection criteria are determined 
in samples of GEANT-4 based [3] Monte Carlo (MC) sim- 
ulated signal decays, where we include the effects of elec- 
tromagnetic radiation from the final-state charged parti- 
cles using the PHOTOS simulation package [l5| . 

We compare the performance of our simulation with 
a scalar QED calculation [ll| resummed to all orders of 
o;qed- Among events selected by the \AE\ < 150 MeV 
requirement, the MC simulation and Ref. predict dif- 
ferent fractions of events with photons with energy be- 
low 2.6 MeV, the soft photon energy cut-off used in our 
simulation (see Tab. We therefore reweight the AE 
distributions for each mode to account for this differ- 
ent fraction of radiating events and use these reweighted 
distributions in the final maximum likelihood fit. The 
difference in event yields obtained with the original dis- 
tributions and with the reweighted ones is used to evalu- 
ate the associated systematic error, and it is found to be 
negligible. 

As explained in Ref. , while taking into account ra- 
diative corrections, one needs to be careful to quote the 
results in such a way that the radiation effects can be 
disentangled. In principle, it would be necessary to se- 
lect B candidates with a specified maximum amount of 
O(100 MeV) photon energy in the final state, a quantity 
that is difficult to reconstruct with the BXBAR detector. 
Instead, we define our data sample by selecting on AE, 
an observable that can be related to the maximum al- 



lowed total energy of the photons, £™ ax . The chosen 
AE window allows for the presence of radiated photons 
with total energy up to 150 MeV ± (AE), where the av- 
erage value of AE, (AE), differs for each mode, due to 
the pion mass hypothesis being assigned to all tracks. As 
the 7r + 7r - events are centered at AE ~ MeV, while the 
K + tt~ and K + K~ distributions are shifted by —45 MeV 
and —91 MeV, respectively, the corresponding energy re- 
quirements on the radiated photons are iJ™ ax = 150, 
105 and 59 MeV for tt + h~ , K + tt~, and K + K~ , respec- 
tively. The smearing of AE due to finite momentum res- 
olution leads to a small difference between the number of 
events that satisfy the AE requirement and the number 
of events that satisfy an equivalent £™ ax requirement. 
We use the MC simulation to evaluate the associated 
systematic error on the branching fractions from this dif- 
ference. 

In addition to signal 7r + 7r~, K + tt~, and (possibly) 
K + K~ events, the selected data sample includes back- 
ground from the process e + e~ — > qq (q — u,d,s,c). Ac- 
cording to the MC simulation, backgrounds from other 
B decays are small relative to the signal yields (< 1%), 
and are treated as a systematic uncertainty. We use an 
unbinned, extended maximum-likelihood (ML) fit to ex- 
tract simultaneously signal and background yields in the 
three topologies (irir, Kit, and KK). The fit uses the 
discriminating variables Toes, AE, the Cherenkov angles 
of the two tracks, and a Fisher discriminant J- , based on 
the momentum flow relative to the h + h'~ thrust axis of 
all tracks and clusters in the event, excluding the h + h'~ 
pair, as described in Ref. 0]. The likelihood for event j is 
obtained by summing the product of the event yield Ni 
and probability Vi over the signal and background hy- 
potheses i. The total likelihood for a sample of N events 
is 



\ i / j L i 



(2) 



The probabilities Vi are evaluated as the product 
of the probability density functions (PDFs) with pa- 
rameters Ui, for each of the independent variables 
xj = {mss, AE,T ,6+ ,6~}, where 0+ and 0~ are the 
Cherenkov angles for the positively- and negatively- 
charged tracks, respectively. We check that the variables 
are almost independent. The largest correlation between 
the Xj is 13% for the pair (toes, AE), and we have con- 
firmed that it has a negligible effect on the fitted yields. 
For both signal and background, the K ±r K T yields are 
parameterized as N K ±^ = Nk-k (1 T Ak-k) /2, and we 
fit directly for the total yield Nkk and the asymmetry 
■Akti- The result for Ak-k is used only as a consistency 
check and does not supersede our previously published 
result 

The eight parameters describing the background 
shapes for rn.Es, AE, and T are allowed to vary freely 
in the ML fit. We use a threshold function [l6| for toes 
(one parameter), a second-order polynomial for AE (two 
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TABLE II: Summary of results from the ML fit for the yields. 
The subscript b refers to background. For the nominal fit, we 
use a double Gaussian for the signal AE PDF, as described in 
the text. We also show, for comparison purposes, the results 
using a single Gaussian, which corresponds to an analysis that 
ignores FSR effects. 



TABLE III: Summary of relative systematic uncertainties on 
signal yields. For the K + K~ yield we show the absolute 
uncertainty. The total uncertainties are calculated as the sum 
in quadrature of the individual contributions. 



Source tt+tt' (%) K + tt~ (%) K+ K~ (events) 



mES 


0.2 


0.4 


1.3 


AE 


0.1 


0.0 


0.3 


signal T 


2.9 


1.5 


2.8 


bkgd T 


0.5 


0.2 


5.9 


9 C quality 


0.2 


0.1 


0.4 


Fit bias 


2.2 


0.9 


1.3 


B bkgd 


0.8 


0.2 


< 0.1 



Total 3.8 1.8 6.8 



Parameter Nominal Fit Ignoring FSR 





485 ± 35 


469 ± 34 


Nktv 


1656 ± 52 


1634 ± 52 


Ak-iv 


-0.136 ±0.030 


-0.135 ±0.030 


N kk 


3± 13 


5± 13 




32983 ± 194 


32998 ± 194 


NbKir 


20778 ± 169 


20801 ± 169 


AbKir 


0.002 ± 0.008 


0.002 ± 0.008 


N bKK 


13358 ± 126 


13356 ± 126 



parameters), and a sum of two Gaussian distributions for 
T (five parameters). For the signal shape in to-es, we use 
a single Gaussian distribution to describe all three chan- 
nels and allow the mean and width to vary in the fit. 
For AE, we use the sum of two Gaussian distributions 
(core ± tail), where the core parameters are common to 
all channels and are allowed to vary freely, and the tail 
parameters are determined separately for each channel 
from the reweighted MC simulation (explained above), 
and fixed in the fit. For the signal shape in T , we use 
an asymmetric Gaussian function with different widths 
below and above the mean. All three parameters are de- 
termined from MC simulation and fixed in the maximum- 
likelihood fit. The 8 C PDFs are obtained from a sample 
of approximately 430000 D*+ D°7r+ (D° R'-tt+) 
decays reconstructed in data, where K~/tt + tracks are 
identified through the charge correlation with the tt + 
from the D* + decay. We construct the PDFs separately 
for K + , K~ , 7r + , and n~ tracks as a function of momen- 
tum and polar angle using the measured and expected 
values of 8 C , and its uncertainty. We use the same PDFs 
for tracks in signal and background events. 

Table [TTI summarizes the fitted signal and background 
yields, and Ktt charge asymmetries. We find a value of 
Aktx consistent with our previously published result 
and a background asymmetry consistent with zero. The 
signal yields are slightly higher than the values reported 
in Ref. [3j due to the removal of the At selection require- 
ment and the addition of the radiative tail in the signal 
AE PDF. In order to quantify the effect of FSR on the 
fitted yields, we perform a second fit using a single Gaus- 
sian for the AE PDF, allowing the mean and width to 
vary. The results are shown in the second column of Ta- 
ble HH where we find that ignoring FSR lowers the tttt 
yield by 3.4% and the Ktt yield by 1.3%. 

As a crosscheck, in Fig.[T]we compare the PDF shapes 
(solid curves) to the data using the event-weighting tech- 
nique described in Ref. (l7| . For each plot, we perform a 
fit excluding the variable being plotted and use the fitted 
yields and covariance matrix to determine the relative 



probability that an event is signal or background. The 
distribution is normalized to the yield for the given com- 
ponent and can be compared directly to the assumed 
PDF shape. We find excellent agreement between the 
data and the PDFs. Figure [2] shows the likelihood ra- 
tio Cs/^Ci for all 69264 events in the fitted sample, 
where Cs is the likelihood for a given signal hypothesis, 
and the summation in the denominator is over all sig- 
nal and background components in the fit. We find good 
agreement between data (points with error bars) and the 
distributions obtained by directly generating events from 
the PDFs (histograms). 

Systematic uncertainties on the branching fractions 
arise from uncertainties on the selection efficiency, sig- 
nal yield, and number of BB events in the sample. Un- 
certainty on the efficiency is dominated by track recon- 
struction efficiency (1.6%) and by the uncertainty on FSR 
(1.3% for tttt, 1.4% for Ktt and 2.9% for KK), which is 
evaluated assuming 100% uncertainty on the smearing 
effect on AE. 

Other systematic uncertainties on selection efficiency 
are those due to requirements on the quality of the 9 C 
measurement (1.0% for tttt, 0.8% for Ktt and 0.5% for 
KK) and on event topology (1.1%). Uncertainty on the 
fitted signal yields is dominated by the shape of the sig- 
nal PDF for T (2.9% for tttt, 1.5% for Ktt) and potential 
bias (2.2% for tttt, 0.9% for Ktt) in the fitting technique, 
as determined from large samples of MC-simulated sig- 
nal events and a large ensemble of pseudo-experiments 
generated from the PDF shapes. Uncertainties due to 
imperfect knowledge of the PDF shapes for toes, AE, 
and 9 C are all less than 1%. Tables [TTT] and [IV] summa- 
rize the uncertainties on the signal yields and branching 
fractions, respectively. 

Tablc[V]summarizes the results for the charge-averaged 
branching fractions. For comparison, we use the efficien- 
cies and signal yields determined under the assumption 
of no FSR and find B{B° -> tt+tt") = 5.0 x 10~ 6 and 
B(B° -> K+tt-) = 18.0 x 10~ 6 , which are consistent 
with our previously published results [7[. We determine 
the upper limit for the signal yield for K + K~ using a 
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FIG. 1: Data distributions (points with error bars) of tubs, AS, and T for signal 7r + 7r~ (a,b,c), signal K + tv~ (d,e,f) and 
background for the three channels (g,h,i), using the weighting technique described in the text. Solid curves represent the 
corresponding PDFs used in the fit. The distribution of AE for signal K + tt~ events is shifted due to the assignment of the 
pion mass for all tracks. 



TABLE IV: Summary of relative systematic uncertainties on 
yields, efficiencies, and number of BB pairs. For the K + K~ 
yield we show the absolute uncertainty. The total uncertain- 
ties are calculated as the sum in quadrature of the individual 
contributions. 



Source 


7T + 7T 


K + TY' 


K+K~ 


yields 


3.8% 


1.8% 


6.8 events 


efficiency 


2.5% 


2.5% 


3.5% 


N BE 


1.1% 


1.1% 


1.1% 


Total 


4.7% 


3.3% 


see text 



Bayesian procedure that assumes a fiat prior on the num- 
ber of events. The upper limit is given by the value of 
N Q for which J No C max dN/ J Q °° C max dN = 0.90, corre- 
sponding to a one-sided 90% confidence interval. Here, 
£ max is the likelihood as a function of the K + K~ yield TV, 
maximized with respect to the remaining fit parameters. 
We find No — 25.4, and the upper limit on the branching 
fraction is calculated by increasing the signal yield up- 
per limit and reducing the efficiency by their respective 
total errors (Table [TV)) . For the purpose of combining 
with measurements by other experiments, we also evalu- 
ate the central value for the branching fraction and find 
B(B° -> K+K-rvy) = (4 ± 15 ± 8) x lO" 8 . 

Although we cannot directly measure the non- 
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FIG. 2: (Color online) Distribution of the likelihood ratio Cs/^d, where Cs is the likelihood for each event to be a signal 
tt + tt~ (left), K + ix~ (middle), or K + K~ (right) event. The points with error bars show the distribution obtained on the fitted 
data sample, while the histograms show the distributions obtained by generating signal (dark shaded, red) and background 
(light shaded, yellow) events directly from the PDFs. 



TABLE V: Summary of branching fraction results. We give 
signal yields Ns, total detection efficiencies (e) and branching 
fractions Be-,, where the subscript _B 7 serves as a reminder of 
the dependence on the cut on soft photon energy as explained 
in the text. The errors are statistical and systematic, respec- 
tively, and the upper limit on B° — > K + K~n^ corresponds 
to the 90% confidence level. 



Mode N s e(%) Be^W^) 

7T+7T" 485 ±35 ±18 41.8 ±0.2 ±1.0 5.1 ± 0.4 ± 0.2 

K + TY~ 1656 ± 52 ± 30 40.5 ± 0.2 ± 1.0 18.1 ±0.6 ±0.6 

K+K~ 3.2 ±12.9 ± 7 39.0 ± 0.3 ± 1.4 < 0.5 (90% C.L.) 



radiative, or "bare" branching fractions, due to the in- 
trinsic and unavoidable features of QED, they can be 
extrapolated from our measurements by employing the- 
oretical calculations, such as those found in Ref. [ill ]. 
The results for these bare branching fractions for the 
three channels are shown in Table IVI1 and the central 
value for the bare K + K~ branching fraction is B°(B° — > 
K + K~) = (4± 15 ±8) x 1CT 8 . We stress the importance 
of being able to disentangle radiation effects from the ex- 
perimental measurements, as a meaningful comparison 
between theory and experiment can be performed only 
in terms of the bare quantities. Likewise, bare quanti- 
ties should be used when combining measurements from 
different experiments. 

In summary, we have presented updated measurements 
of charge-averaged branching fractions for the decays 
B" -> 7T+7T- and B° K+tt', with FSR effects taken 
into account. We find that the branching fractions are 
a few percent higher when the effect of FSR is included 
in the calculation of the efficiency and signal yield de- 
termination. This difference should be taken into ac- 
count when comparing with previous measurements of 



these quantities 0, 0, 0, EH that do not include these 
effects. In order to perform the most meaningful com- 

TABLE VI: Summary of experimental branching fractions, 
Be-! , with a defined cut on soft photon energy, together with 
the electromagnetic correction factor G(i?™ ax ) and the eval- 
uated "bare" branching fractions (non radiative), B° . The 
errors on branching fractions are statistical and systematic 
respectively; the error on G(i?™ ax ) is taken as the difference 
between its value at fj, = M n and fi = M p . 



Mode Be-,(W~ u ) G(£™ ax ) g u (10~ b ) 

tt^tP 5.1 ±0.4 ±0.2 0.937 ±0.005 5.5 ± 0.4 ± 0.3 

K + TT~ 18.1 ±0.6 ±0.6 0.947 ±0.005 19.1 ±0.6 ±0.6 

K+ K ~ < 0.5 (90% C.L.) 0.952 ± 0.005 < 0.5 (90% C.L.) 



parison, we also evaluated the bare branching fractions 
for the three channels, as explained in Ref. 11]. Our 
results are consistent with current theoretical estimates 
from different models Q. We find no evidence for the 
decay B° — > K + K~ and set an upper limit of 5.0 x 1CP 7 
at the 90% confidence level. 
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